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ABSTRACT 

Monoallelic expression within a gene family is 
found in pathogens exhibiting antigenic variation 
and in mammalian olfactory neurons. Trypanosoma 
brucei, a lethal parasite living in the human blood- 
stream, expresses variant surface glycoprotein 
(VSG) from 1 of 15 bloodstream expression sites 
(BESs) by virtue of a multifunctional RNA polymerase 
I. The active BES is transcribed in an extranucleolar 
compartment termed the expression site body (ESB), 
whereas silent BESs, located elsewhere within the 
nucleus, are repressed epigenetically. The regulatory 
mechanisms, however, are poorly understood. Here 
we show that two essential subunits of the basal 
class I transcription factor A (CITFA) predominantly 
occupied the promoter of the active BES relative to 
that of a silent BES, a phenotype that was maintained 
after switching BESs in situ. In these experiments, 
high promoter occupancy of CITFA was coupled 
to high levels of both promoter-proximal RNA 
abundance and RNA polymerase I occupancy. 
Accordingly, fluorescently tagged CITFA-7 was 
concentrated in the nucleolus and the ESB. Because 
a ChlP-seq analysis found that along the entire BES, 
CITFA-7 is specifically enriched only at the promoter, 
our data strongly indicate that monoallelic BES tran- 
scription is activated by a mechanism that functions 
at the level of transcription initiation. 

INTRODUCTION 

Monoallelic expression of a large gene family is a phenom- 
enon found in parasites that exhibit antigenic variation of 



cell surface proteins to evade the mammalian host's 
immune system (1). It also occurs in mammahan olfactory 
sensory neurons, which express a single member of the 
large odorant receptor (OR) gene family (2). While 
monoallelic gene expression of a gene family is typically 
regulated at the level of transcription, the mechanisms 
enabhng transcription of a single gene/allele are not well 
understood. 

Trypanosoma brucei is a vector-borne lethal parasite in 
sub-Saharan Africa that lives freely in the human blood- 
stream by virtue of a dense coat of variant surface glyco- 
protein (VSG) expressed from a single gene. The uniform 
coat shields invariant membrane proteins from immune 
recognition (3) and is periodically altered by switching 
VSG expression to a different gene. While the parasite 
harbors more than 2000 VSG genes and pseudogenes, 
VSG expression is restricted to 1 of 15 telomeric blood- 
stream expression sites (BESs) (4). A BES is a tandem 
array of several expression site-associated genes 
(ESAGs) followed by a single VSG gene, all of which 
are polycistronically transcribed by a multifunctional 
RNA polymerase (pol) I from a promoter that is located 
40-60 kb upstream of the VSG gene (5). The active BES 
is transcribed outside of the nucleolus (6) in the 
extranucleolar expression site body (ESB), a DNase- 
resistant compartment that appears to limit productive 
transcription to a single BES (7). 

Silencing of gene family members is mediated by epi- 
genetic mechanisms. Silent OR genes in mice are marked 
with histone methylations H3K9me3 and H4K20nie3 that 
are characteristic of constitutive heterochromatin (8). 
Plasmodium falciparum, the parasite causing cerebral 
malaria, exhibits monoallelic expression of the var gene 
family, which encodes variants of the P. falciparum 
erythrocyte membrane protein 1 (9). This protein is 
expressed on the surface of an infected erythrocyte and 
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mediates cytoadherence to endothelial host receptors 
keeping the infected cell in the peripheral vascular 
system. Deficiency of the histone deacetylase Sir2 in 
P. falciparum led to activation of silent var genes and 
loss of monoallelic var gene expression (10,11). 
Accordingly, knockdown of T. brucei RAPl, encoding a 
telomeric protein known to induce epigenetic gene 
silencing, led to derepression of aU known BESs, detection 
of additional extranucleolar RNA pol I foci and forma- 
tion of mixed VSG ceU surface coats (12). In addition, 
several other factors have been implicated in BBS silencing 
because their gene knockdowns resulted in derepression 
of silent BESs; these comprise the methyltransferase 
DOTIB (13), the chromatin remodelers TbISWI (14) 
and FACT (15), the histone deacetylase DAC3 (16), the 
nucleoplasmin-like protein NLP (17), histone HI (18), 
the origin recognition complex (19), histone H3 and the 
histone chaperones CAP- lb and ASFIA (20) as well 
as the mini-chromosome maintenance-binding protein 
MCM-BP (21). Moreover, nucleosomes are enriched in 
the chromatin of silent BESs but are depleted in that of 
the active BES (22,23). 

While these studies demonstrated the importance of epi- 
genetic factors in restricting VSG expression to a single 
gene, the mechanism of transcriptional silencing is yet to 
be determined. Importantly, in the cited studies, mature 
VSG mRNA levels from derepressed BESs remained ~ 10- 
10000-fold lower than that of the active VSG. BES dere- 
pression by RAPl and DOTIB silencing did not affect the 
high expression level of the active VSG (12,13), indicating 
that derepressed BESs cannot effectively compete with 
VSG expression from the active BES. On the other 
hand, silencing either SPT16, encoding a FACT subunit, 
or NLP did reduce the abundance of active VSG tran- 
script (15,17). However, as shown in the NLP study, 
promoter-proximal gene expression of the active site was 
much less affected, suggesting again that in the dere- 
pressed state, BES promoters cannot efficiently compete 
for transcription factors from the active promoter. 

We have recently characterized the multi-subunit class I 
transcription factor A (CITE A) as an essential factor 
for RNA pol I transcription of the active BES, the large 
ribosomal gene unit (RRNA) and the procychn gene 
units, which encode the major cell surface proteins of 
insect-stage, procyclic form (PF) trypanosomes (24). 
Characterization of CITFA revealed seven apparently 
trypanosomatid-specific proteins, termed CITFA- 1 to -7, 
and the dynein light chain DYNLLl (24,25). The purified 
complex bound to the 67-bp-long BES promoter in a gel 
shift assay requiring the integrity of both known promoter 
elements for efficient binding. Moreover, depletion of 
CITFA-2 or CITFA-7 from extract or inactivation of 
CITFA-2 by an anti-CITFA-2 antibody abohshed RNA 
pol I transcription in vitro from BES, rRNA gene (RRNA) 
and procychn gene promoters, as determined by 121-146- 
bp-long primer extension products. This identified CITFA 
as a basal promoter-binding transcription initiation factor 
for RNA pol I in trypanosomes (24,25). Accordingly, 
silencing of CITFA-2 and of CITFA-7 resulted in rapid 
loss of rRNA and mRNA of the active VSG gene in 
bloodstream forms (BFs) and in trypanosome death in 



culture. Here, anti-CITFA-2 and anti-CITFA-7 chroma- 
tin inimunoprecipitation (ChIP) assays showed that 
CITFA predominantly occupied the active BES 
promoter relative to that of a silent BES, a phenotype 
maintained after switching BESs in situ. High CITFA oc- 
cupancy at the active BES promoter correlated with high 
promoter occupancy by RNA pol I and with a high 
promoter-proximal RNA level. Accordingly, as shown 
for the RRNA loci, CITFA-7 silencing considerably 
reduced both promoter-proximal RNA pol I occupancy 
and RNA abundance. Together these data strongly 
indicate that monoallelic VSG transcription is regulated 
at the level of transcription initiation. Moreover, 
fluorescently tagged CITFA-7 was found to be 
concentrated in the nucleolus and the ESB supporting 
the notion that CITFA has a role in BES activation. 
Because a ChlP-seq analysis revealed that among BES 
sequences high CITFA-7 occupancy is restricted to the 
promoter, our results raise the possibihty that the 
sequestration of CITFA limits productive transcription 
initiation by RNA pol I to the nucleolus and the ESB. 

MATERIALS AND METHODS 

DNAs 

pCITFA-7-PTP-NEO and the RNAi stem-loop vector 
pCITFA-7-T7-stl were described previously (25). pCITFA- 
7-eYFP-NEO was generated by inserting 353 bp of the 
CITFA-7 coding region (position 56 to position 408) into 
the Apal and NotI sites of pC-eYFP-NEO (26). pCITFA- 
7-PTP-BLA was made by replacing the neomycin 
phosphotransferase with the blasticidin S deaminase ORE, 
pBLA-PTP-CITFA-2 and pBLA-PTP-RPB6z were likewise 
generated from pPURO-PTP-CITFA-2 (24) and pPURO- 
PTP-RPB6Z (27). pPURO-mCherry-RPB6z was obtained 
by replacing the sequence of the composite PTP tag 
(protein C epitope, tobacco etch virus [TEV] protease site, 
tandem protein A domain) with the mCherry fruit fluores- 
cent protein coding region. Allele deletions were achieved by 
transfecting polymerase chain reaction (PCR) products in 
which 100 bp of 5' and 3' gene flanks were fused to the 
hygromycin phosphotransferase coding region. DNA oligo- 
nucleotides that were used in competitive and quantitative 
PCR (qPCR) are listed in Supplementary Table SI. 
Semiquantitative reverse transcriptase-PCR (RT-PCR) was 
performed using cycle numbers which were empirically 
determined to be within the linear amplification range for 
each oligonucleotide pair. To demonstrate that consensus 
oligonucleotides co-amplify promoter sequences from differ- 
ent BESs and RRNA repeats, we sequenced PCR amphfica- 
tion products of genomic DNA and detected 11 out of 14 
and 2 out of 3 known single-nucleotide polymorpliisms of 
BES (4) and RRNA promoters, respectively (Supplementary 
Figure SI and data not shown). 

Cells 

Culturing and stable transfection of BFs and PFs were 
carried out as previously described (25,28). BF cells were 
grown with either G418, puromycin, hygromycin, 
phleomycin or blasticidin at final concentrations of 2.5, 
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0.2, 1, 1 and 2|ig/ml, respectively, unless specified other- 
wise. Double-strand CITFA-7 RNA expression was 
induced with 2|.ig/ml doxycycline. For in situ BBS 
switching, 1.8 x 10 Np/nP cells were washed three times 
with prewarmed medium, resuspended in medium to 
2x10 cells/ml, divided into six ahquots and each 
aliquot grown in antibiotic-free medium to a 13-ml 
culture of 1 x 10'' cells/ml. After dilution with an equal 
volume of medium, cells were plated in 1 ml aliquots of 
either lO^g/ml puromycin or 100|ig/ml G418. Stable 
clonal cell fines were maintained at eitfier 1 |.ig/ml puro- 
mycin or 10 |.ig/ml G418. Assays of drug sensitivity of ceUs 
following BES switchings were as described in (13), 
namely, 10 ml of 1 x 10 cells/ml were grown in either 
1 ^g/ml puromycin or lOO^g/ml G418 for 48 h, growth 
was then monitored for 72 h. 

Fluorescence microscopy 

Microscope slides and coverslips (2.5x2.5 cm) were 
washed in methanol in a heated ultrasonic water bath 
for 20min, rinsed thoroughly with deionized water and 
washed with 100% ethanol for 20min. Following a final 
rinse in water, the shdes and coverslips were air-dried. 
A square chamber was made on each sfide using strips 
of labeling tape. Cells in exponential growth were 
centrifuged for 4min at 800^ and room temperature, 
washed in cold medium and resuspended in cold BBSG 
buffer (50 mM Bicine, 77 mM glucose, 50 mM NaCl and 
5mM KCl, pH 8.0) to a density of 1 x 10^ ceUs/nl. CeUs 
were fixed by adding 4% paraformaldehyde to a final con- 
centration of 0.05% for 1 min after which they were 
pelleted at 800^ for 45 s. Following two washes with 
BBSG buffer, the cell pellet was resuspended in BBSG 
buffer to a final concentration of 4 x 10 cells/|il and 4,6- 
diamidino-2-phenyhndol (DAPI) was added to a final con- 
centration of 0.04|ig/ml. Twenty-five microliters of cell 
suspension were spread on the chamber and covered 
with 125 |il of 1% low gelling agarose kept warm at 
37°C (Sigma, St. Louis, MO, USA), and a coverslip was 
placed on top to close the chamber. Slides were stored at 
4°C for ~16h for complete cell fixation and pictures were 
acquired within 24 h from the time of embedding, using an 
Axiovert 200 microscope (Zeiss, Thornwood, NY, USA). 
Images were processed with Axiovision software (Zeiss) 
and ImageJ. Whole cell and ESB enhanced yellow fluor- 
escent protein (eYFP) signals were quantified using 
ImageJ and corrected by subtracting the mean pixel inten- 
sity of the surrounding region of each cell analyzed. 
Individual cell analysis was restricted to cells with only 
one nucleus and one kinetoplast. 

Chip and ChlP-seq analyses 

Chip experiments were carried out as previously described 
(29), resulting in average chromatin fragments of 200- 
400 bp. Analysis of Np/nP cell fines required larger chro- 
matin fragments to accommodate amplification of the se- 
lectable marker regions downstream of the transcription 
initiation site (TIS). DNA fragments of ~600bp were 
obtained with medium sonication power using a 
Bioruptor UCD-200 (Diagenode, New York, NY, USA) 



and a reduced sonication time from 25 to 8 min (30 s on/ 
30 s off). PTP-tagged CITFA-7, CITFA-2 and RPB6z 
were immunoprecipitated with an anti-ProtA rabbit poly- 
clonal serum (Sigma) and untagged RPA31 with a rat 
anti-RPA31 immune serum (27). Negative control precipi- 
tations were conducted with either rabbit or rat 
nonspecific immune sera. The percent IP was calculated 
relative to the input material and corrected by subtracting 
the percent IP of the negative control assays. The fold 
enrichment was obtained by dividing the amount of 
precipitated DNA by the amount of DNA from control 
precipitations. To generate the library for ChlP-seq 
analysis, 1 x lO' cells were processed to produce 10 ng of 
precipitated DNA. In this case, chromatin was sonicated 
for 100 min (ten 10 min of 30 s on/30 s off with 5 min 
coofing on ice in between) on high setting to obtain 
smafier fragments ranging from 100 to 250 bp in size. 
qPCR quality controls and setup were analyzed as 
described previously (29). For the preparation of ChlP- 
seq libraries, 10 ng of input DNA and 10 ng of DNA 
precipitated with CIFTA-7 were end-repaired by 
incubating the DNA with 1 x T4 figase buffer (NEB), 
0.4 mM dNTPs, 0.5^1 of End Repair Enzyme Mix 
(NEB, 30 \il of 3 U/nl T4 DNA polymerase, 6 ^1 of 1 U/ 
|.il Klenow Fragment, 30|il of lOU/^l T4 DNA 
Polynucleotide Kinase) in a total volume of 50|il for 
30 min at 20°C. The end-repaired DNA fragments were 
cleaned up using AmPure XP beads (Beckman Coulter) 
as suggested by Illumina. To generate a protruding 3' A 
base for adapter ligation, the end-repaired DNA was 
mixed with 2^1 of 10 x NEB buffer 2, of 4mM 
dATP and 0.5 |il of Klenow 3' to 5' exo minus enzyme 
(5U/|.il, NEB) and incubated for 30 min at 37°C. Next, 
the DNA was figated to 1 |il of forked DNA adapters 
(0.2 |iM) using the Quick Ligation Kit (NEB) according 
to the manufacturer's specifications. To prepare the 
forked adapters, forward (AATGATACGGCGACCAC 
CGAGATCTACACTCTTTCCCTACACGACGCTCTT 
CCGATC*T) and reverse DNA oligonucleotides (/5Phos/ 
GATCGGAAGAGCACACGTCTGAACTCCAGTCA 
CGCCAATAT CTCGTATGCCGTCTTCTGCTT*G or/ 
5Phos/GATCGGAAGAGCACACGTCTGAACTCCA 
GTCACCTT GTAATCTCGTATGCCGTCTTCTGCTT 
*G) were suspended at 100 ^iM in lOmM Tris-HCl, pH 
7.8, 0.1 mM EDTA, pH 8.0, 50mM NaCl, mixed at equal 
concentrations, heated to 95°C and annealed by slowly 
coofing the mix to room temperature (asterisk denotes a 
phosphorothioate bond). Adaptor-annealed DNA was 
cleaned up using AmPure XP beads and amplified using 
the KAPA HIFI HotStart Ready Mix (Kapabiosystems) 
and PCR primers 5'-AATGATACGGCGACCACC 
GA*G-3' (forward primer) and 5'-CAAGCAGAAGAC 
GGCATACGA*G-3' (reverse primer) for five cycles as 
suggested by Kapabiosystems. The PCR product was 
gel-purified and, again, PCR-aniplified (15 cycles). The 
final PCR product was cleaned up using the MinElute 
PCR Purification Kit (Qiagen) and gel purified on a 1% 
agarose gel followed by a cleanup using the MinElute Gel 
extraction Kit (Qiagen). The quality and concentration of 
the final libraries was determined using a BioAnalyzer 
2100 (Agilent). 
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Libraries from input and immunoprecipitated DNA were 
sequenced using a HiSeq2000 sequencer (llluniina). For tlie 
input DNA, we obtained 75 167 347 sequence reads of 
whicli 73.6% could be mapped (45.6% mapped to exactly 
1 site). For the immunoprecipitated DNA, we obtained 
16722060 sequence reads of which 59.8% could be 
mapped (39.8% mapped to exactly one site). The 
sequence reads (100 bp) were annotated based on the 
T. brucei 427 genome (version 4.1) to which we added 
the sequences of the 14 BESs available at http://www. 
sanger.ac.uk/resources/downloads/protozoa/trypanosoma- 
brucei.html. Standard and fast gapped read alignments 
were performed using Bowtie version 2 (30). To normaUze 
the read coverage between the libraries, the read coverage 
was divided by the total number of aligned reads in each 
library and multiplied by the total number of aligned reads 
of the smaller library. 

RNA analysis 

Total RNA was prepared by the hot phenol method 
as specified previously (27). For pre-mRNA analysis, 
total RNA was reverse transcribed with random 
hexanucleotides (Roche) as primers and Superscript II 
reverse transcriptase (Invitrogen) following the manufac- 
turer's recommendations, and cDNA was amplified with 
at least one oligonucleotide hybridizing outside the 
sequence of the mature niRNA. 

RESULTS 

CITFA-7 occupancy of BES and RRNA promoters 

To further investigate the transcriptional roles of the 
CITFA complex in RNA pol I transcription, we 
analyzed CITFA-7 occupancy of RRNA and BES pro- 
moters in BFs. Previously, we have generated the BE 
cell line TbC7ee that expresses CITFA-7 with a C- 
terminal FTP tag and no untagged CITFA-7 (Figure 1) 
(25). The large composite PTP tag which consists of a 
protein C epitope, a TEV protease site and dual protein 
A (ProtA) domains (31) did not affect the functionahty of 
CITFA-7 because CITFA-7 is essential for trypanosome 
viabihty and TbC7ee cells prohferated as efficiently as 
wild-type cells (25). Using a ChlP-grade, polyclonal anti- 
ProtA antibody (32) and a nonspecific rabbit immune 
serum for normalization, we precipitated chromatin 
from TbC7ee cells. According to qPCR analysis using 
consensus ohgonucleotides for BES and for RRNA pro- 
moters, both promoter sequences were clearly enriched in 
the precipitate in comparison with the P— a tubulin 
intergenic region, which served as a negative control 
(Figure 1). DNA fragments within the ESAG7 and 18S 
rRNA coding regions, which extended to ~1.6 and 
~1.2kb downstream of the BES and RRNA TISs, respect- 
ively, precipitated much less efficiently than their pro- 
moters. Because we have previously found that RNA 
pol I cross-hnked equally to the RRNA promoter and to 
the 18S rRNA coding region (29), these results strongly 
indicated that CITFA-7 does not associate with 
elongating RNA pol I and that it binds only to the 
promoter region, thus substantiating the role of CITFA 
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Figure 1. CITFA-7 occupies specifically RNA pol I promoters. Top 
panel, schematic drawing of the CITFA-7 locus in cell line TbC7ee, 
which exclusively expresses CITFA-7 with a C-terminal PTP fusion. 
The open reading frames of CITFA-7, of //yC" and NEO^, and the 
PTP tag are indicated by an open box, striped boxes and a black box, 
respectively. Smaller gray boxes represent gene flanks for RNA pro- 
cessing signals. The PTP tag sequence was fused to an endogenous 
CITFA-7 allele by targeted integration of pCITFA-7-PTP as indicated. 
Bottom panel, anti-ProtA ChIP with TbC7ee cells expressing CITFA-7- 
PTP. Precipitated DNA was analyzed by qPCR of the consensus BES 
promoter (BES prom), a conserved region of the ESAG7 coding region 
{ESAGl ORF), the consensus RRNA promoter {RRNA prom), 18S 
rRNA coding region (18S rDNA) and the P — a tubuhn intergenic 
region (P — a tub). The precipitation percentage relative to the input 
material was corrected by subtracting corresponding values from 
negative control precipitations that were carried out with a comparable 
nonspecific antiserum. 



as a transcription initiation factor. Most interestingly, 
RRNA promoter DNA was consistently ~ 5-fold more 
enriched than BES promoter DNA. Because the trypano- 
some harbors multiple copies of these promoters, this 
result suggested that CITFA-7 occupied relatively more 
RRNA than BES promoters and, therefore, raised the pos- 
sibility that among BES promoters it primarily bound the 
promoter of the active BES. 

CITFA-7 predominantly occupies the active BES 
promoter relative to a silent BES promoter 

The previously described BF cell fine PN13 in which the 
silent BESl and the active BES 17 were marked with a 
puromycin N-acetyltransferase gene (PURO^) and a 
neomycin phosphotransferase (NEO^) gene, respectively 
(13), was used to directly test the binding of CITFA-7 to 
promoters of active and silent BESs. Both genes were 
introduced at corresponding positions 240 bp downstream 
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of the TIS, locations that allowed for the selection of in 
situ switches between these two BESs and for specific 
DNA amplification downstream of the BESl and BES17 
promoters. To enable anti-CITFA-7 ChIP experiments in 
these cells, we stably integrated plasmid CITFA-7-PTP- 
BLA into one CITFA-7 allele and termed the correspond- 
ing line C7-Npl, indicating that these trypanosomes are 
resistant to neomycin (upper case N), sensitive to puro- 
mycin (lower case p) and represent the basal fine from 
which nP2 and Np3 cell lines were directly derived by con- 
secutive in situ BES switching experiments (Figure 2A). 

Anti-ProtA IP of chromatin prepared from C7-Npl 
cells was first analyzed by a competitive PCR assay in 
which a common sense primer was paired with two anti- 
sense ohgonucleotides that were complementary to either 
the NEO^ or the PURO^ coding region (Figure 2A). 
While the NEO^ band was only slightly stronger than 
the PURO^ band in DNA isolated from total chromatin 
or from the control IP, anti-ClTFA-7 precipitated chro- 
matin was strongly enriched in NEO^ DNA derived from 
the active BES 17 (Figure 2B). Incubation of C7-Npl cells 
without antibiotics and subsequent selection with puro- 
mycin combined with limiting dilution, revealed independ- 
ent cell fines that had switched to the expression of BES 1 . 
A corresponding ChIP analysis of one of these fines 
(C7-nP2) showed that now the BESl -derived PURO^ 
band was highly enriched over the NEO^ band, 
demonstrating that high CITFA-7 occupancy correlated 
with BES activation (Figure 2B). 

To quantitatively assess CITFA-7 occupancy at BESl/ 
BES 17 promoters in these cell lines, we carried out qPCR 
analysis of precipitated DNA. To ensure specific amplifi- 
cation of BESl and BES 17 sequences, we chose ofigo- 
nucleotide pairs that hybridized to the PURO^ and 
NEO^ coding regions, respectively. As shown for the 
original C7-Npl cell line and six independently derived 
C7-nP2 clones and C7-Np3 clones, ClTFA-7 occupied 
predominantly the active BES promoter relative to that 
of the silenced BES (Figure 2C). This pattern was not 
due to continuous selective antibiotic pressure because it 
remained unchanged when the antibiotics were removed 
for ten generations before the ChIP experiment. We did, 
however, consistently detect above background levels of 
CITFA-7 at the promoter of the silenced BES, indicating 
that the transcription factor does interact with silent BES 
promoters to some extent. 

Although we specifically amplified the active BES in 
these experiments, the precipitation percentage did not 
increase over that of the assays of the consensus BES 
promoter (compare Figure 2C with Figure 1). This has 
two fikely explanations: firstly, TbC7ee cells exclusively 
expressed tagged CITFA-7, whereas C7-Np/nP lines co- 
expressed wfid-type and tagged CITFA-7. Secondly, in the 
first experiment (Figure 1), we amplified the BES 
promoter, which is the CITFA binding site, whereas the 
specific NEO^ and PURO'^ amplification products of the 
subsequent experiment (Figure 2A) extended to positions 
482 and 545 downstream of the TIS, respectively, thus 
amplifying only a fraction of the precipitated DNA. 
Accordingly, when we calculated the fold enrichment 
over negative control precipitations, the active BES 17 in 



C7-Npl cells exhibited a similarly high enrichment as the 
consensus RRNA promoter, a nearly 8-fold higher enrich- 
ment than that of the inactive BESl promoter, and an 
approximately four times higher enrichment than that of 
the consensus BES promoter in TbC7ee cells, confirming 
high CITFA-7 occupancy of the active BES promoter 
(Figure 2D). 

BES activation is regulated at the level of 
transcription initiation 

Our previous in vitro transcription results strongly 
indicated that CITFA is a basal, DNA-binding transcrip- 
tion initiation factor that mediates RNA pol I recruitment 
to its cognate promoters. The purified CITFA complex 
required both BES promoter elements for efficient DNA 
binding and depletion of CITFA from extract resulted in 
the loss of RRNA and BES promoter transcripts that were 
only 127 and 146 nucleotides long, respectively (24;25). 
Hence, high CITFA occupancy should lead to high 
promoter-proximal transcription rates and to high RNA 
pol I occupancy at the active BES promoter. To assess 
promoter-proximal transcription levels at the active 
versus the silenced BES, we quantified the abundances 
of pre-TVE'O^ and pie-PURO'^ mRNA relative to tubufin 
pre-mRNA in C7-Npl cells as weU as in one C7-nP2 and 
one C7-Np3 cell fine (Figure 2E). Because in these cell 
lines NEO'^ and PURO^ RNA processing is directed by 
aldolase gene flanks (13), we paired a sense oligonucleo- 
tide that targets the aldolase 5' gene flank upstream of 
the known aldolase splice site with either a NEO^- or 
PURO^-specific antisense ofigonucleotide to amplify 
random hexamer-derived cDNA. As expected, pre- 
mRNA of the selectable marker from the active site was 
16 (Npl), 40 (nP2) and 24 (Np3) times more abundant 
than tubulin pre-mRNA and exceeded the selectable 
marker pre-mRNA from the inactive BES 59- to 79-fold 
(Figure 2E). 

Because this finding suggested a higher transcription 
initiation rate from the active BES promoter than from 
the inactive BES promoter, we wanted to determine the 
promoter-proximal occupancy of RNA pol I at these sites. 
We generated analogous to the C7-Npl cell fine 
(Figure 2 A), an R6z-Npl cell line in which we introduced 
the functional N-terminal PTP fusion (27) to the RNA pol 
I-specific subunit RPB6z in PN13 cells. The corresponding 
Chip experiment demonstrated that RNA pol I occu- 
pancy was ~ 17-fold higher at the active BES 17 
promoter than at the silent BESl promoter (Figure 2F). 

Finally, we confirmed that CITEA-7 silencing is func- 
tionally linked to the loss of promoter-proximal RNA pol 
I occupancy and the reduced abundance of promoter- 
proximal pre-RNA. Because BESs in our extant single 
marker cells for gene knockdowns were unmarked, we 
were unable to specifically investigate the active BES. 
Instead, we focused our analysis on the RRNA promoter 
region because CITFA-7 is equally important for RRNA 
and BES transcription (25), and it is fikely that most RRNA 
repeats are active in rapidly dividing, cultured trypano- 
somes. For these ChIP experiments we used a polyclonal 
antibody directed against the essential RNA pol I-specific 



Nucleic Acids Research, 2014, Vol. 42, No. 5 3169 



BES17 



BES1 



cell line C7-Np1 

P'^^ NEO" ESAG7 



■*-2 

PURCf 



ESAG7 



CITFA-7 



CITFA-7-PTP BLA" 




— NEO" 

— PURO" 

— NEO" 
— PURO" 



0.1- 



o 0.05- 
o 



0 



0.25TT- 



■2 0.2-H 

Q. 
O 
OJ 

CL 

0.15-H 

0) 

o 
a3 

■a 

(U 

-•— ' 

o 



□ BESMNECf 
■ BES1 PURO" 

□ p-a tub 



C7-Np1 -nP2.1 -nP2.2 -nP2.3 -Np3.1 -Np3.2 -Np3.3 



F RNA pol I occupancy 



60 



_ 45 
I 30 

(D 

o 



15 



C7-Np1 



RRNA 
prom 



BES 
prom 



BES17 BES1 

A/EO" PURO" 



0) 100 
o 

•a 



< 
z 
ir 
E 

Q. 

> 



^o■■■ 



1... 



■ BES1 7 NEO" 

■ BES1 PURCf 



I 

I 



C7-Np1 -nP2 -Np3 



o 



0.3 



g. 
o 

(i> 
i_ 

Q. 
+-< 

0) 

o 

i_ 

a> 

Q- 

■o 

0) 

o 

o 
o 



0) 0.2 



0.1 



*** 
I 1 



□ BES17 A/EO" 
■ BES1 PURCf 
n P-a tub 



cell line RBz-Npl 
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PCR amplifying NEO'^- and Pf/^O^-specific sequences in total DNA (Inp), control IP (ctrl) and anti-ProtA IP (IP). (C) qPCR analysis of anti- 
CITFA-7-PTP Chip assays, carried out with clonal cell lines obtained from in situ switching events between BES 17 and BESl, using NEO"- or 
/"(/T^O^-specific primer pairs. The difference between BES17 and BESl occupancies were statistically analyzed by paired two-tailed student's /-tests. 
One, two and three asterisks indicate /'-values that are <0.05, 0.01 and 0.001, respectively. (D) Calculation of the fold enrichment over negative 
control precipitations for the ChIP experiments shown in the bottom panel of Figure 1 and the C7-Npl cell line shown in Figure 2C (BES17 NEO^ 
and BESl PURO"). (E) Reverse transcription-qPCR analysis of NEO'^ or PURO" pre-mRNA abundances in C7-Np/nP cell lines. The values for 
these relative abundances were standardized with the tubulin pre-mRNA level, which was set to 1. (F) qPCR analysis of anti-PTP-RPB6z ChIP with 
R6z-Npl cells. Results shown in (E) and (F) were statistically analyzed as described above. 



subunit RPA31 (27). This antibody effectively precipitated 
the RRNA promoter and the 18S rRNA coding region in 
noninduced cells, while the (5— a tubulin intergenic region, 
which is transcribed by RNA pol II, was not precipitated 
above the control precipitation with a comparable nonspe- 
cific polyclonal antibody (Figure 3A). On silencing CITFA- 
7 for 42 h, RNA pol I occupancy decreased, on average, by 
83% at the RRNA promoter and by 84% in the 18S rRNA 
coding region. Similarly, promoter-proximal 18S pre- 
rRNA abundances were reduced, on average, by 91.5% 
in C/Ti^^- 7-silenced cells as analyzed by qPCR of 



random hexamer-derived cDNA. The corresponding 
semiquantitative qPCR analysis shows the drop of pre- 
rRNA abundance and the concomitant increase of 
tubulin pre-mRNA (Figure 3B). The latter finding is con- 
sistent with our previous observation that the massive loss 
of rRNA in C/7Ty4-7-silenced cells led to an increase of all 
nonaffected RNAs in total RNA preparations of equal 
amounts (25). These results confirmed the pivotal role of 
CITFA in RNA pol I transcription initiation. 

Taken together, we have shown here that high CITFA-7 
occupancy of the active BES promoter coincides with high 
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Figure 3. CITFA-7 silencing affects RRNA promoter-proximal RNA pol I occupancy and pre-18S rRNA abundance. (A) Anti-RPA31 ChIP assays 
were conducted with noninduced bloodstream trypanosomes (BFs) and with BFs in which CITFA-7 was silenced for 42 h. Precipitated chromatin was 
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levels of promoter-proximal RNA from this BES and high 
RNA pol I occupancy in the promoter region. In addition, 
we demonstrated for the RRNA gene array that depletion 
of ClTFA-7 strongly reduces RNA pol I occupancy in 
the promoter region and causes significant decreases of pro- 
moter-proximal RNA abundance. Therefore, we concluded 
that monoallelic BES transcription in T. brucei is regulated 
at the level of transcription initiation. 

BES promoter occupancy by CITFA-2 parallels that 
of CITFA-7 

The high occupancy of the active BES promoter by 
ClTFA-7 suggested that the whole CITFA complex is 
bound there. To confirm this notion, we repeated the 
Chip analysis with the essential CITFA-2 subunit. 
Similar to the CITFA-7 analysis, we generated cell line 
C2-Npl by tagging CITFA-2 N-terminally with PTP in 
PN13 cells. We previously showed that PTP-CITFA-2 is 
fully functional (24). The competitive PCR analysis 
demonstrated even more strikingly than in the CITFA-7 
analysis the predominant association of CITFA-2 with the 
active BES promoter (Figure 4A). qPCR analyses of anti- 
PTP-CITFA-2 Chip assays conducted in C2-Npl cells 
and in independently derived, switched C2-nP2 and C2- 
Np3 hues showed that CITFA-2 was even more distinctly 
associated with the active BES than CITFA-7 (Figure 4B). 
While in case of CITFA-7, on average 3.47 times more 
DNA was precipitated from the active versus the silent 
BES promoter, this fold difference in CITFA-2 occupancy 
was, at 6.24, significantly higher (P = 2.4 x 10~^, two- 
tailed, unpaired student's ^-test assuming unequal 
variance). Previously, by sedimenting tandem affinity- 
purified CITFA complex through a linear sucrose 
gradient, we obtained evidence that trypanosomes 
possess two CITFA complexes: a more abundant 
complex that lacks CITFA-2 and DYNLLl and a 
shghtly faster sedimenting full complex (25). Because 
CITFA-2 is indispensable for CITFA function in RNA 
pol I transcription from BES, RRNA and GREET 
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Figure 4. High CITFA-2 occupancy of the active BES promoter. 

(A) Anti-PTP-CITFA-2 ChIP assays with C2-Npl and C2-nP2 cell 
lines analyzed by competitive PCR of NEO" and PURO" sequences. 

(B) Corresponding qPCR analysis of C2-Np/pN cell lines. Statistical 
analysis was by paired two-tailed Student's r-tests. 

procyclin gene promoters, it appears that binding of 
CITFA-2 to the complex activates CITFA. Therefore, 
the more distinct association of CITFA-2 with the active 
BES promoter may be directly related to the fact that 
precipitating CITFA-2 specifically pulled down the 
active CITFA complex. Together, these data strongly 
indicated that, among BES promoters, CITFA is predom- 
inantly bound to the active BES promoter. 
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CITFA-7 is concentrated in the nucleolus and the ESB 

The efficiency of the RNA pol I transcription system is 
due to its high initiation rate, which, as was shown for the 
yeast ribosomal gene unit (33), can lead to an RNA pol I 
density of one actively transcribing enzyme complex for 
every ~100 bp. Accordingly, detection of RNA pol I in the 
ESB has been straightforward [reviewed in (34)]. On the 
other hand, it was demonstrated in the human system that 
the basal, promoter-binding RNA pol I transcription ini- 
tiation factors SLl and UBF can support multiple rounds 
of transcription in vitro without dissociating from the 
promoter (35), which suggests that relatively few tran- 
scription factors are required to ensure maximal RNA 
pol I transcription. Because the ESB appears to harbor a 
single BES with either one or two BES promoters (4), it 
was not certain whether CITFA would be detectable in 
this compartment. However, in a previous study, when 
we used indirect immunofluorescence microscopy to co- 
localize the PTP-tagged RNA pol I subunit RPB6z and 
HA-tagged CITFA-7, we occasionally detected a CITFA- 
7 signal outside the nucleolus that co-localized with 
RPB6z, suggesting that the ESB contains detectable 
amounts of the transcription factor (25). To enhance the 
sensitivity of the RPB6z/CITFA-7 co-locahzation, we 
created a BF cell line in which, by targeted plasmid inte- 
grations into the RPB6z and CITFA-7 loci, the mCherry 
protein and the eYFP were N- and C-terminally fused 
to RPB6z and CITFA-7, respectively (Figure 5A). 
Fluorescence microscopy showed that in 39% of cells 
(« = 155) mCherry-RPB6z was detected only in the nucle- 
olus, which is apparent in DAPI stains as a spherical 
structure of low DNA density (Figure 5B, long arrow). 
In 41% of cells, mCherry-RPB6z localized to the nucle- 
olus as well as to a single extranucleolar spot that was 
clearly separate from the nucleolus and, therefore, repre- 
sented the ESB. Like the nucleolus, the ESB consistently 
exhibited a reduced DAPI stain (Figure 5B, short 
arrow). In the remaining 20% of cells, two extranucleolar 
mCherry-RPB6z signals were detectable. These numbers 
were in accordance with previous RNA pol I-based ESB 
identifications (7,13). Although the nucleolar staining 
patterns of CITFA-7-eYFP and mCherry-RPB6z were 
not completely congruent, which could be expected from 
a transcription initiation factor and an RNA pol, these 
proteins exhibited extensive co-locahzation in the nucle- 
olus of each cell analyzed. Most importantly, in 91% of 
the cells with one or two extranucleolar spots, we could 
detect co-localization of CITFA-7-eYFP in afl these 
spots (Figures 5B, C, two additional examples are shown 
in Supplementary Figure S2A). Although after our 
previous study we had anticipated that CITFA-7 would 
be detectable in the ESB, the rather robust signal of 
CITFA-7-eYFP in this compartment was nevertheless 
surprising because it suggested a much larger accumula- 
tion of the transcription factor in the ESB than is required 
for RNA pol I transcription initiation from a single BES. 
Quantification of eYFP fluorescence revealed that the ESB 
signal contributed 4.5% to that of a whole cell (standard 
deviation of 1.3%, n = 10). Because an immunoblot 
comparison of recombinant CITFA-7 with endogenous 
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Figure 5. Co-localization of mCherry-RPB6z and CITFA-7-eYFP. (A) 
Schematic of the RPB6~ and CITFA-7 loci showing the targeted 
pPURO-mCherry-RPB6z and pCITFA-7-eYFP-NEO integrations in 
the BF cell line used here. (B) Images from top to bottom show 
phase contrast, DAPI stain (bar: 10|im; n, nucleus; k, kinetoplast; 
long and short arrows indicate nucleolus and ESB, respectively), 
mCherry-RPB6z, CITFA-7-eYFP (pseudo-colored in green) and an 
image merge. Images on the left are unmodified, whereas red and 
green signals in the enlarged images (panels in the middle and on the 
right) were recolored and enhanced in brightness and contrast. (C) 
Statistical analysis of 155 individual cells. Cells were categorized by 
the number of extranucleolar RNA pol I (mCherry-RPB6z) spots. 
Cells in which mCherry-RPB6z and CITFA-7-eYFP exhibited clear 
co-localization in the nucleolus and all detectable extranucleolar spots 
were scored as cells with complete co-localization. 



CITFA-7 of whole cell lysates uncovered that a BF 
trypanosome harbors ~4000 CITFA-7 molecules 
(Supplementary Figure S2B), we calculated that the ESB 
contains ~180 CITFA-7 molecules. Together, these data 
show that CITFA-7 is concentrated in both the nucleolus 
and the ESB, which strongly supports our finding of 
predominant CITFA occupancy of the active BES 
promoter. 
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Peak(s) of CITFA-7 occupancy within a BES is 
restricted to its promoter region(s) 

Our findings of high CITFA occupancy of the active BES 
promoter and the presence of ~200 CITFA-7 molecules in 
the ESB raised the question of how this factor is 
concentrated in this compartment. We have previously 
shown in vitro that the CITFA complex is a sequence- 
specific DNA binding factor that requires both BES 
promoter elements for efficient binding (24). However, it 
was possible that BESs harbor multiple CITFA binding 
sites or that CITFA-7 remains bound to RNA pol I during 
transcription elongation, which would lead to an enrich- 
ment of CITFA-7 along the entire BES. To determine the 
presence of CITFA along the entire BES, we conducted 
a genome-wide CITFA-7-PTP ChlP-seq analysis and 
aligned the reads to the genome of T. brucei brucei All 
including the 14 BES sequences, which were available 
for this strain (4). To account for sequencing biases and 
alignment artifacts, the number of sequencing reads 
obtained from the IP was normalized with the number 
of sequencing reads obtained from the input material. 
Reads (100 nt) were aligned to the most similar genomic 
sequence or, if the target site was not unique, randomly 
distributed among all possible target sites. Thus, the 
analysis was limited by two constraints: firstly, BES 
promoter regions are nearly identical, differing only in a 
few single nucleotide polymorphisms, which made unam- 
biguous alignment of most reads to particular BES pro- 
moters impossible. Hence, reads from the active site were 
presumably distributed among most BESs. Secondly, for 
our analysis, we used cells that originated from the Tc221 
clone of the 427 strain (36) that had switched expression 
from the original VSG221 gene to an unknown VSG 
gene (kind gift from Dr Peter Overath, MPI Tiibingen, 
Germany). As we found out during our analysis, the 
BES promoter regions of this strain differed in a few pos- 
itions from the pubhshed 427 strain (4), preventing the 
unambiguous identification of individual BESs according 
to the known promoter sequences (data not shown). 
Nevertheless, sequence reads could be aUgned along all 
BESs, indicating that the Tc221 BES sequences still 
resembled the published ones enough to enable a mean- 
ingful analysis. 

Figure 6 shows the analysis of BES 15 (also annotated as 
telol26), which has two promoters of highly similar 
sequence (PI and P2). After normalizing for the different 
number of reads obtained with input and immunopre- 
cipitated DNA, peaks of reads were observed in several 
BES regions, most of them in the 70-bp repeat regions, 
the VSG gene and the VSG pseudogene (Figure 6, second 
panel). However, these peaks were also observed in input 
DNA that was not enriched by IP (Figure 6, third panel), 
which suggests that they were formed owing to an 
underrepresentation of these repetitive sequences in the pub- 
lished genome. Two peaks in BES 15, however, were absent 
in the input. They precisely mapped to the two BES 15 pro- 
moters and, therefore, are the result of ChIP enrichment. 
Accordingly, when we calculated the fold enrichment of 
Chip reads over input reads, 13.6- and 9.2-fold enrichments 
were obtained for PI and P2, respectively. Closer inspection 



of the 67-bp-long P2 promoter showed that the highest 
number of reads and the highest folds of enrichment pre- 
cisely mapped to all seven windows representing the 
promoter sequence (Figure 6, zoomed-in panels). The 
number of reads on either side of the promoter rapidly 
declined, suggesting that the transcription factor occupies 
only the promoter and does not follow the polymerase 
after transcription initiation. This notion was backed up 
by our finding that RPAl, the largest subunit of RNA pol 
I, did not co-immunoprecipitate with CITFA-7-PTP in a 
DNA-independent manner (Supplementary Figure S3) and 
by the fact that subunits of RNA pol I or of CITFA never 
cross-purified in several tandem affinity purifications 
(24,25,27,37,38). Hence, we concluded that, within a BES, 
CITFA-7 occupancy is restricted to the promoter region 
and that there are no additional binding sites that could 
explain CITFA concentration in the ESB. This result is in 
close agreement with our standard ChIP results, which 
showed that CITFA-7 occupancy was strongly decreased 
downstream of the BES promoter (Figure 1). Given that 
the single 67-bp BES promoter is unlikely to support 
binding of multiple CITFA complexes, and that CITFA is 
not found away from this promoter, our data raise the inter- 
esting possibility that CITFA is concentrated in the nucle- 
olus and the ESB independent of DNA and RNA pol I 
association, thereby restricting high RNA pol I transcrip- 
tion initiation rates to these two compartments. 

CITFA-7 was enriched to some extent on aU BES pro- 
moters analyzed (Supplementary Figure S4 and data not 
shown). This hkely resulted from reads of the active BES 
being distributed among other BESs. In support of this 
notion, we found more reads at RRNA promoters than on 
BES promoters. Within the nine annotated RRNA pro- 
moters, the peak window accumulated an average of 
1233 normalized reads, which is 4.25-fold more than the 
BES 15 promoter peak, suggesting that reads from the 
active BES were aligned to silent BESs owing to 
sequence identities. Moreover, we observed several-fold 
differences in the number of normalized reads assigned 
to the different BES promoters that were phylogenetically 
Hnked (Supplementary Figure S5). The combined reads 
for the seven promoter windows ranged from 1 54 for the 
single BESl promoter to 2148 for the single BES5 
promoter. Promoters whose sequence closely resembled 
those of the similar BES15/BES5 promoters obtained 
high numbers of reads, whereas more distantly related 
BES promoters obtained much fewer reads. These 
findings are in agreement with CITFA being predomin- 
antly bound to the promoter of the active BES. 



DISCUSSION 

VSG expression in T. brucei is an ideal system to study the 
mechanism of monoallelic gene expression because RNA 
pol I transcription of the active BES in the ESB is spatially 
separated from RRNA transcription in the nucleolus. We 
took advantage of this fact and showed that CITFA-7, an 
essential subunit of the CITFA complex, is concentrated 
in the nucleolus and the ESB, substantially co-localizing 
with RNA pol I in these compartments. Consistent with 
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Figure 6. CITFA-7 enrichment at BES15 {telol26). For BES15, the fold enrichment of DNA immunoprecipitated with CITFA-7 compared with 
input DNA (first panel) and the norinalized read coverages obtained from the CITFA-7 IP (second panel) and the input material (third panel) are 
shown. The read coverage was normalized to account for differences in library size (see 'Materials and Methods' section) and calculated for a 
window size of 11 bp. Orange boxes represent ESAGs identified by numbers (e.g. seven stands for ESAG7) and of VSG genes (* identifies 
pseudogenes), green boxes indicate promoters, and black boxes depict 70-bp repeats. For the zoomed-in panels, sequences for each window are 
specified. The two BES promoter elements and the TIS position are identified by red lettering and the remaining promoter sequence by blue lettering. 
The bars corresponding to the promoter sequence are depicted in blue. The green bar at the bottom represents the promoter region annotation in the 
genome data bank. 
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this observation, we found that CITFA-7 and the previ- 
ously characterized, essential subunit CITFA-2 were pre- 
dominantly associated with the promoter of the marked 
active BBS versus that of a marked silent BES. In situ 
switching between these two BBSs reversed the promoter 
occupancy pattern of the two CITFA subunits, which es- 
tablished that BBS activation is coupled with high CITFA 
interaction. Because CITFA is a promoter-binding tran- 
scription initiation factor that is absolutely required for 
RNA pol I transcription in vivo and in vitro (24,25), high 
CITFA occupancy should lead to high rates of RNA pol I 
recruitment and promoter-proximal transcription rates. 
Accordingly, we found the RNA pol I-specific subunit 
RPB6z to be highly abundant in the promoter region of 
the active BBS, and promoter-proximal RNA levels of the 
active BES to exceed those of the silent BBS ~70-fold. 
Furthermore, knockdown of ClTFA-7 led to a sharp 
decrease of RNA pol I occupancy at the RRNA 
promoter and the downstream 18S RRNA coding 
region, as well as a high reduction of the promoter- 
proximal pre-rRNA. Together, these data show that 
monoallelic VSG expression is controlled, at least par- 
tially, at the level of transcription initiation and indicate 
that this control is mediated by binding of a functional 
CITFA complex to the active BBS promoter. This is the 
first demonstration in any system of antigenic variation 
that basal transcription initiation factor occupancy 
differs between the active and the silent genes. 

Although CITFA predominantly occupied the active 
BBS promoter, we clearly observed above background 
CITFA occupancy at the silent BBS promoter. One pos- 
sibihty is that the marked silent BES in our study was not 
completely silenced but remained in a 'pre-active' state. 
Previously, as a model to explain co-expression of two 
BBSs under antibiotic selection, it was proposed that 
between silent and active BBSs there is a pre-active state 
that precedes BES activation (39). The pre-active state of a 
BBS is characterized by close spatial proximity to the 
active BBS and thought to enable rapid back and forth 
switching between these two BBSs under selective 
pressure. Because we switched repeatedly between BESl 
and BES 17, it is possible that the silenced BES was locked 
in the pre-active state giving it greater access to CITFA 
than a completely silenced BES. On the other hand, our 
detection of CITFA at silent BES promoters is consistent 
with previous studies demonstrating that a low level of 
transcription does initiate at silent BBSs but does not 
reach the promoter-distal VSG gene (40,41). A previous 
expression analysis of the promoter-proximal ESAG6 gene 
revealed that ~20% of all ESAG6 mRNA were derived 
from silent BESs, whereas 80% originated from the active 
site (40). Taking 15 BESs into account, this implies an 
~ 56-fold stronger transcription initiation rate from the 
active BBS than from any of the silent sites. This fold 
difference closely reflects our pre-mRNA-based estimate 
of an ~70-fold higher initiation rate at the active BBS 
promoter. However, cloning and sequencing of RT-PCR 
products along the complete BBS revealed a more diverse 
picture and promoter-proximal expression from silent 
BBSs varied strongly between different T. brucei strains 
(41). Nevertheless, in all cases, promoter-proximal RNA 



levels derived from the active BBS did exceed those from 
any silent BES manifold. Together, these results confirm 
that activation of a BBS not only requires ablation of 
telomeric silencing but also a mechanism that leads to 
higher transcription initiation rates. 

Knowledge about proteins that preferentially interact 
with the promoter regions of silenced or active BBSs is 
limited. While we do know that silent BESs are occupied 
by nucleosomes that are absent from the active BES 
(22,23), it is not known whether this difference is a conse- 
quence of transcriptional activity or represents an actual 
control mechanism (42). In agreement with the latter, a 
recent study showed that the HMG (high mobility group) 
box protein TDPl, which is thought to be a general archi- 
tectural chromatin protein, exhibited high occupancy of 
the active BBS and RRNA gene units, thereby likely 
facihtating their transcription (43). Additionally, histone 
HI (18) and NLP (17) were found to be mainly associated 
with silent BBS promoters, which suggest that they 
have specific roles in reducing transcription from these 
promoters. However, it remains unclear whether the dif- 
ferential distributions of these proteins among BESs are 
regulatory in nature or whether they reinforce the 
activated or silenced status of a BES. Nevertheless, these 
studies support our model that there are control mechan- 
isms that operate at BES promoters. 

Our data show that the active BBS promoter recruits 
CITFA more efficiently than a promoter of a silent BBS. 
One possible explanation for this difference is that repres- 
sive chromatin interferes with CITFA binding to silent 
BBS promoters. This notion is supported by experiments 
that showed that silencing of DAC3, ISWl and NLP led to 
an upregulation in promoter-proximal gene expression at 
silent BESs (14,16,17). However, in these cases the dere- 
pressed expression level was not compared with that of the 
active BES and it is hkely that it remained considerably 
lower. The high transcription level of the active BES (44) 
hkely impedes upregulation of all silent BBSs to the same 
level. If BBS activation is achieved by derepression only, 
one would expect that general BES derepression would 
effectively impede transcription of the active site, owing 
to competition for CITFA and RNA pol I. Thus far, 
promoter-proximal transcription was only assessed in 
the NLP knockdown, which led to a moderate ~ 3 -fold 
reduction of mRNA from a selectable marker gene that 
was introduced downstream of the active BES promoter 
(17). Because our results suggest that the rate of RNA pol 
I transcription initiation is limited by CITFA availability, 
it will be important to analyze whether BES derepression 
leads to higher CITFA occupancy. 

An alternative explanation for the predominant associ- 
ation of CITFA with the active BES would be a mechan- 
ism that confines this basal transcription factor to the 
nucleolus and the BSB, limiting its availability and high 
RNA pol I transcription initiation rates to these compart- 
ments. Because we did not find evidence for a direct inter- 
action between CITFA-7 and RNA pol I and CITFA 
appeared to be cross-hnked only to the promoter region 
in our ChlP-seq analysis, the high concentration of 
CITFA-7 in these compartments seemed to be independ- 
ent of RNA pol I and the capability of CITFA to bind 
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DNA. A model of confined CITFA availability would be 
in agreement with previous analyses of VSG expression. 
The spatial proximity of co-expressed BESs observed by 
Chaves et al. (39) could reflect the need to acquire CITFA 
from the ESB for selectable marker expression. This 
model would also be consistent with previous findings 
showing that functional replacement of a EES promoter 
with an RRNA promoter left BES regulation intact (45) 
because CITFA also enables RRNA promoter transcrip- 
tion (24). Finally, it concurs with results obtained on 
silencing of RAPl, which led to a significant increase of 
detectable RNA pol I accumulations outside the nucleolus 
but did not alter the transcription efficiency of the 
promoter-proximal BES regions (12). These findings 
imply that transcription elongation of RNA pol I along 
BESs is mainly blocked epigenetically and that, in accord- 
ance with our results, VSG activation occurs at the BES 
promoter region. Lastly, these results show that RNA pol 
I accumulation at derepressed BESs is not sufficient to 
drive promoter-proximal BES transcription, suggesting 
that RNA pol I sequestration is not regulating BES tran- 
scription, whereas confinement of CITFA may do so. 

It should be noted here that CITFA subunits are 
conserved among all trypanosomatids including those that 
do not have a multifunctional RNA pol I. Except for some 
structural differences between RRNA and BES promoters 
(46), RRNA transcription initiation in the nucleolus and 
BES transcription in the ESB appear to function equiva- 
lently. Nevertheless, because protein sequestration to the 
nucleolus is a well-established phenomenon (47), it is 
possible that the ESB is a nucleolus-like structure that 
offers interaction sites for the CITFA complex. As discussed 
previously, the dynein component of CITFA is a likely can- 
didate to facilitate factor sequestration (24). However, 
DYNLLl binding motifs, which are conserved from yeast 
to mammals, were not detected in CITFA sequences and, 
thus, will have to be determined experimentally. 

Spatial confinement of an essential gene expression 
factor may have relevance beyond T. brucei's RNA pol I 
system. In P. falciparum, the active var gene is transcribed 
by RNA pol II in a perinuclear site and activation of 
a silent var gene requires relocation from the hetero- 
chromatic nuclear periphery to this euchromatic site 
possibly to interact with a sequestered var gene-specific 
transcription factor [reviewed in (9)]. Similarly, DNA 
FISH analysis in the mouse system showed that the H 
enhancer element co-localized with OR genes on dif- 
ferent chromosomes, suggesting that nuclear positioning 
of the active OR gene is important in this system as 
well (48). 

In summary, our results favor a bipartite model of 
monoallelic VSG expression in which BES repression is 
achieved epigenetically and the expressed BES is not only 
derepressed by chromatin changes but also activated at the 
promoter, a process that entails recruitment of CITFA. 



SUPPLEMENTARY DATA 

Supplementary Data are available at NAR Online. 



ACKNOWLEDGEMENTS 

We are grateful to Drs George Cross and Luisa 
Figueiredo (Rockefeller University) for providing the 
PN13 cell line and the BES switching protocol, to 
Dr Ziyin Li (University of Texas, Houston) for construct 
pC-eYFP-NEO, to Drs Ji Yu and Dongmyung Oh 
(UConn Health Center) for help with fluorescence micros- 
copy, to Ms Bao Nguyen for help with plasmid construc- 
tion and to Dr Justin Radolf and Justin K. Kirkham for 
critical reading of the manuscript. 



FUNDING 

This work was funded by the National Institutes of Health 
grant [ROl AI059377 to A.G.]; the Young Investigator 
Program of the Research Center of Infectious Diseases 
(ZINF) of the University of Wiirzburg, Germany, and 
the Human Frontier Science Program (to T.N.S.). 
Funding for open access charge: National Institutes of 
Health [ROl AI059377 to A.G.]. 

Conflict of interest statement. None declared. 



REFERENCES 

1. Horn.D. and McCuUoch.R. (2010) Molecular mechanisms 
underlying the control of antigenic variation in African 
trypanosomes. Curr. Opin. Microbiol., 13, 700-705. 

2. Magklara,A. and Lomvardas,S. (2013) Stochastic gene expression 
in mammals: lessons from olfaction. Trends Cell Biol., 23, 
449-456. 

3. Schwede,A-, Jones, N., Engstler,M. and Carrington.M. (2011) The 
VSG C-terminal domain is inaccessible to antibodies on live 
trypanosomes. Mol. Biochem. ParasiloL, 175, 201-204. 

4. Hertz-Fowler,C., Figueiredo,L.M., Quail,M.A., Becker.M., 
Jack:son,A., Bason,N., Brooks,K., Churcher,C., Fahkro.S., 
Goodhead,!. et al. (2008) Telomeric expression sites are highly 
conserved in Trypano.soma hriicei. PLoS One, 3, e3527. 

5. GiinzkA., Bruderer,T., Laufer,G., Scliimanski,B., Tu.L.C, 
Chung,H.M., Lee,P.T. and Lee,M.G. (2003) RNA polymerase I 
transcribes procyclin genes and variant surface glycoprotein gene 
expression sites in Trvpano.sonm hrucei. Eukarvot. Cell, 2, 
542-551. 

6. Chaves,!., Zomerdijk.J., Dirks, M. A., Dirks, R.W., Raap,A.K. and 
Borst,P. (1998) Subnuclear localization of the active variant 
surface glycoprotein gene expression site in Trypano.soma hrucei. 
Proc. Natl Acad. Sci. USA, 95, 12328-12333. ' 

7. Navarro,M. and Gull,K. (2001) A pol I transcriptional body 
associated with VSG mono-allelic expression in Trypanosoma 
hrucei. Nature, 414, 759-763. 

8. Magklara,A., Yen.A., Colquitt,B.M., Clowney,E.J., Allen,W., 
Markenscoff-Papadimitriou,E., Evans.Z.A., Kheradpour,?., 
Mountoufaris.G., Carey.C. et al. (2011) An epigenetic signature 
for monoallelic olfactory receptor expression. Cell, 145, 555-570. 

9. ScherfA., Lopez-Rubio,J.J. and Riviere,L. (2008) Antigenic 
variation in Plasmodium falciparum. Anmi. Rev. MicrohioL, 62, 
445-470. 

10. Duraisingh.M.T., Voss,T.S., Marty.A.J., Duffy,M.F., Good,R.T., 
Thompson.J.K., Freitas-Junior,L.H., Sclierf A., Crabb,B.S. and 
Cowman,A.F. (2005) Heterochromatin silencing and locus 
repositioning linked to regulation of virulence genes in 
Plasmodium falciparuin. Cell, 121, 13-24. 

11. Tonkin,C.J., Carret.C.K., Duraisingh.M.T., Voss.T.S., Ralph,S.A., 
HommefM., Duffy,M.F., Silva,L.M., Scherf A., Ivens.A. et al. 
(2009) Sir2 paralogues cooperate to regulate virulence genes and 
antigenic variation in Plasmodium falciparum. PLoS Biol., 1, e84. 



3176 Nucleic Acids Research, 2014, Vol. 42, No. 5 



12. Yang,X., Figueiredo,L.M., Espinal,A., Okubo,E. and Li,B. (2009) 
RAPl is essential for silencing telomeric variant surface 
glycoprotein genes in Trypanosoma bnicei. Cell, 137, 99-109. 

13. Figueiredo,L.M., Janzen,C.J. and Cross,G.A. (2008) A histone 
methyltransferase modulates antigenic variation in African 
trypanosomes. PLoS Biol., 6, el61. 

14. Hughes, K., Wand,M., Foulston,L., Young,R., Harley.K., 
Terry,S., Ersfeld,K. and Rudenko,G. (2007) A novel ISWI is 
involved in VSG expression site downregulation in African 
trypanosomes. EMBO J., 26, 2400-2410. 

15. Denninger,V., Fullbrook,A., Bessat,M., Ersfeld,K. and 
Rudenko,G. (2010) The FACT subunit TbSptl6 is involved in 
cell cycle specific control of VSG expression sites in Trypaiio.ioma 
hrucei. Mol. Microbiol., 78, 459^74. 

16. Wang,Q.P., Kawahara,T. and Horn,D. (2010) Histone 
deacetylases play distinct roles in telomeric VSG expression site 
silencing in African trypanosomes. Mol. Microbiol., 11, 
1237-1245. 

17. Narayanan, M.S., Kushwaha,M., Ersfeld,K., Fullbrook,A., 
Stanne,T.M. and Rudenko,G. (2011) NLP is a novel transcription 
regulator involved in VSG expression site control in Trypano.sonia 
hrucei. Nucleic Acids Res., 39, 2018-2031. 

18. Povelones,M.L., Gluenz,E., Dembek,M., Gull,K. and Rudenko,G. 
(2012) Histone HI plays a role in heterochromatin formation and 
VSG expression site silencing in Trypanosoma brucei. PLoS 
Pathog., 8, el003010. 

19. Tiengwe,C., Marcello,L., Farr,H., Dickens,N., Kelly,S., 
Swiderski,M., Vaughan,D., Gull,K., Barry,J.D., BellS.D. et al. 
(2012) Genome-wide analysis reveals extensive functional 
interaction between DNA replication initiation and transcription 
in the genome of Trypanosoma brucei. Cell Rep., 2, 185-197. 

20. Alsford.S. and Horn,D. (2012) Cell-cycle-regulated control of 
VSG expression site silencing by histones and histone chaperones 
ASFIA and CAF-lb in Trypanosoma brucei. Nucleic Acids Res., 
40, 10150-10160. 

21. Kim,H.S., Park,S.H., Gunzl,A. and Cross,G.A. (2013) MCM-BP 
is required for repression of life-cycle specific genes transcribed by 
RNA polymerase I in the mammalian infectious form of 
Trypanosoma brucei. PLoS One, 8, e57001. 

22. Figueiredo,L.M. and Cross, G.A. (2010) Nucleosomes are depleted 
at the VSG expression site transcribed by RNA polymerase I in 
African trypanosomes. Eukaryot. Cell, 9, 148-154. 

23. Stamie,T.M. and Rudenko,G. (2010) Active VSG expression sites 
in Trypanosoma brucei are depleted of nucleosomes. Eukaryot. 
Cell, 9, 136-147. 

24. Brandenburg,!., Schimanski,B., Nogoceke,E., Nguyen,T.N., 
Padovan,J.C., Chait,B.T., Cross,G.A. and Gunzl,A. (2007) 
Multifunctional class I transcription in Trypanosoma brucei 
depends on a novel protein complex. EMBO J., 26, 4856^866. 

25. Nguyen,T.N., Nguyen,B.N., Lee,J.H., Panigrahi,A.K. and 
GiinzkA. (2012) Characterization of a novel class I transcription 
factor A (CITFA) subunit that is indispensable for transcription 
by the multifunctional RNA polymerase I of Trypanosoma brucei. 
Eukaryot. Cell, 11, 1573-1581. 

26. Li,Z., Umeyama,T. and Wang,C.C. (2009) The Aurora Kinase in 
Trypanosoma brucei plays distinctive roles in metaphase-anaphase 
transition and cytokinetic initiation. PLoS Pathog., 5, el000575. 

27. Nguyen,T.N., Schimanski,B. and Gunzl,A. (2007) Active RNA 
polymerase I of Trypanosoma brucei harbors a novel subunit 
essential for transcription. Mol. Cell. Biol., 27, 6254-6263. 

28. Laufer,G., Schaaf,G., Bollg6nn,S. and Gunzl.A. (1999) In vitro 
analysis of alpha-amanitin-resistant transcription from the rRNA, 
procyclic acidic repetitive protein, and variant surface 
glycoprotein gene promoters in Trypcmosoma brucei. Mol. Cell. 
Biol, 19, 5466-5473. 

29. Park,S.H., Nguyen,T.N., Kirkham,J.K., Lee,J.H. and Gunzl,A. 
(2011) Transcription by the multifunctional RNA polymerase I in 
Trypanosoma hrucei functions independently of RPB7. Mol. 
Biochem. ParasitoL, 180, 35-42. 



30. Langmead,B. and Salzberg,S.L. (2012) Fast gapped-read 
alignment with Bowtie 2. Nat. Methods, 9, 357-359. 

31. Schimanski,B., Nguyen,T.N. and Gunzl,A. (2005) Highly efficient 
tandem affinity purification of trypanosome protein complexes 
based on a novel epitope combination. Eukaryotic Cell, 4, 
1942-1950. 

32. Lee,J.H., Cai,G., Panigrahi,A.K., Dunham-Ems, S., Nguyen,T.N., 
Radolf,J.D., Asturias,F.J. and Gunzl,A. (2010) A TFIIH- 
associated mediator head is a basal factor of small nuclear spliced 
leader RNA gene transcription in early-diverged trypanosomes. 
Mol. Cell. Biol, 30, 5502-5513. 

33. French,S.L., Osheim,Y.N., Cioci,F., Nomura,M. and Beyer,A.L. 
(2003) In exponentially growing Saccharomyces cerevisiae cells, 
rRNA synthesis is determined by the summed RNA polymerase I 
loading rate rather than by the number of active genes. Mol. Cell 
Biol., 23, 1558-1568. 

34. Navarro,M., Penate,X. and Landeira,D. (2007) Nuclear 
architecture underlying gene expression in Trypcmosoma brucei. 
Trends Microbiol, 15, 263-270. 

35. Panov,K.I., Friedrich,J.K. and Zomerdijk,J.C. (2001) A step 
subsequent to preinitiation complex assembly at the ribosomal 
RNA gene promoter is rate limiting for human RNA 
polymerase I-dependent transcription. Mol. Cell Biol., 21, 
2641-2649. 

36. Hirumi,H., Hirumi,K., Doyle,J.J. and Cross,G.A. (1980) In vitro 
cloning of animal-infective bloodstream forms of Trypanosoma 
brucei. Parasitology, 80, 371-382. 

37. Nguyen,T.N., Schimanski,B., Zahn,A., Klumpp,B. and GiinzLA. 
(2006) Purification of an eight subunit RNA polymerase I 
complex in Trypanosoma brucei. Mol. Biochem. ParasitoL, 149, 
27-37. 

38. Walgraffe,D., Devaux,S., Lecordier,L., Dierick,J.F., Dieu,M., Van 
Den,A.J., Pays,E. and Vanhamme,L. (2005) Characterization of 
subunits of the RNA polymerase I complex in Trypanosoma 
brucei. Mol. Biochem. ParasitoL, 139, 249-260. 

39. Chaves,!, Rudenko,G., Dirks,M.A., Cross,M. and Borst,P. (1999) 
Control of variant surface glycoprotein gene-expression sites in 
Trypanosoma brucei. EMBO J., 18, 4846-4855. 

40. Ansorge,I., Steverding,D., Melville,S., Hartmann,C. and 
Clayton,C. (1999) Transcription of 'inactive' expression sites in 
African trypanosomes leads to expression of multiple transferrin 
receptor RNAs in bloodstream forms. Mol. Biochem. ParasitoL, 
101, 81-94. 

41. Vanhamme,L., Poelvoorde,P., Pays,A., Tebabi,P., Xong,H.V. and 
Pays,E. (2000) Differential RNA elongation controls the variant 
surface glycoprotein gene expression sites of Trypcmosoma brucei. 
Mol. Microbiol., 36, 328-340. 

42. Rudenko,G. (2010) Epigenetics and transcriptional control in 
African trypanosomes. Essays Biochem., 48, 201-219. 

43. Narayanan,M.S. and Rudenko,G. (2013) TDPl is an HMG 
chromatin protein facilitating RNA polymerase I transcription in 
African trypanosomes. Nucleic Acids Res., 41, 2981-2992. 

44. Ehlers,B., Czichos,J. and Overath,P. (1987) RNA turnover in 
Trypanosoma hrucei. Mol. Cell Biol., 1, 1242-1249. 

45. Rudenko,G., Blundell,P.A., Dirks Mulder,A., Kieft,R. and 
Borst,P. (1995) A ribosomal DNA promoter replacing the 
promoter of a telomeric VSG gene expression site can be 
efficiently switched on and off in T. hrucei. Cell, 83, 547-553. 

46. Gunzl,A., Vanhamme,L. and Myler,P.J. (2007) Transcription in 
trypanosomes: a different means to the end. In: Barry,J.D., 
Mottram,J.C., McCulloch,R. and Acosta-Serrano,A. (eds), 
Trypanosomes — After the Genome. Horizon Press, Pittsburgh, 
pp. 177-208. 

47. Emmott,E. and Hiscox,J.A. (2009) Nucleolar targeting: the hub 
of the matter. EMBO Rep., 10, 231-238. 

48. Lomvardas,S., Barnea,G., Pisapia,D.J., Mendelsohn,M., 
Kirkland,J. and Axel,R. (2006) Interchromosomal interactions 
and olfactory receptor choice. Cell, 126, 403^13. 



